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Abstract: We calculate the cross section for s-wave neutralino annihilation to 00 
three-body final states below the W~^W~ and ti thresholds. Such three-body chan- '^"^ — ^ 
nels may dominate the annihilation cross section if the neutralino mass is not too 
much less than rrit and mw respectively. Furthermore, because neutrinos produced in 

these channels are much more energetic than those from the bb or r^r^ channels, they | ^ 

can dominate the energetic-neutrino fiuxes from neutralino annihilation in the Sun or 
Earth far below these thresholds and significantly enhance the neutrino signal in certain 
regions of the supersymmetric parameter space. 
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1. Introduction 

It has long been well established that the observed luminous matter in Galactic halos 
cannot account for their total mass. Determination of the identity of this unseen dark 
matter has become one of the most important problems in modern cosmology. Perhaps 
the most promising dark-matter candidate is the neutralino x 5 ^ linear combination 
of the supersymmetric partners of the Z, 7, and Higgs bosons. 

The existence of neutralinos in our halo could be inferred by observation of ener- 
getic neutrinos from annihilation of neutralinos that have accumulated at the core of 
the Sun and/or Earth in detectors such as AMANDA, super-Kamiokande, MACRO, 
and HANUL 0. These neutrinos are produced by decays of neutralino annihilation 
products such as r leptons, c, b, and t quarks, and gauge and Higgs bosons if the 
neutralino is heavy enough. In all cases where the signal is expected to be observable 
by these or subsequent-generation detectors, accretion of neutralinos from the Galactic 
halo onto the Sun or Earth comes into equilibrium with their depletion through annihi- 
lation 0. Therefore, the annihilation rate depends on the rate of capture of neutralinos 
from the halo. 

Although the total annihilation cross section is not needed for flux predictions, the 
branching ratios into the various annihilation products are: The rate for observation 
of neutrinos from decays of various annihilation products may differ considerably. En- 
ergetic neutrinos are much more easier to detect than low-energy ones. For example, 
the flux of energetic neutrinos from decays of b quarks is much smaller than that from 
gauge bosons or top quarks with the same injection energy. 



The best technique for inferring the existence of these neutrinos is to observe an 
upward muon produced by a charged-current interaction in the rock below the detector. 
The rate for observation of energetic neutrinos is proportional to the second moment 
of the neutrino energy spectrum, so it is this neutrino energy moment weighted by the 
corresponding branching ratio that determines the detection rate. 

By now, the cross sections for annihilation have been calculated for all two-body 
final states that arise at tree level (fermion-antifermion and gauge- and Higgs-boson 
pairs). The precise branching ratios depend on numerous coupling constants and su- 
perpartner masses. However, roughly speaking, among the two-body channels the bb 
and T'^T~ final states usually dominate for < mw Neutralinos that are mostly 
higgsino annihilate primarily to gauge bosons if > mw, because there is no s-wave 
suppression mechanism for this channel. Neutralinos that are mostly gaugino continue 
to annihilate primarily to bb pairs until the neutralino mass exceeds the top-quark mass, 
after which the tt final state dominates, as the cross section for annihilation to fermions 
is proportional to the square of the fermion mass. 

Three-body final states arise only at higher order in perturbation theory and are 
therefore usually negligible. However, as we already noted, some two-body channels 
easily dominates the cross section when they are open because of their large couplings, 
for example the W~^W~ for the higgsinos and ti for gauginos. This suggests that their 
corresponding three-body final states can be important just below these thresholds. 

More importantly, as mentioned above, the rate for indirect detection is also pro- 
portional to the second moment of the energy of the muon neutrino. The neutrinos 
produced in these three-body final states are generally much more energetic than those 
produced in b and r decays, so even when the cross sections to these channels are small 
compared with others, they could dominate the indirect-detection rate. 

In this paper, we calculate the cross section for the processes XX ~^ W'^W~* 
Wff' and XX ~^ tt* — > tW~b, and their charge conjugates (heretofore referred as tWb 
and WW* states) in the Vj-e\ limit, where the star denotes virtual particles. Our 
calculation for the WW* is applicable for generic neutralinos, but it is mostly significant 
for neutralinos that are primarily made of higgsino. This is because the gauginos have 
small couplings to W^W~ pairs. As for the tWb final state, our calculation is only 
applicable to the gaugino. If the neutralino is primarily higgsino, then annihilation 
to this final state could also proceed through a WW* intermediate state, which we 
have not included in our calculation. However, if the neutralino is primarily higgsino 
with mw < m^ < mt, then it will annihilate primarily to pairs below the tt 

threshold. Furthermore, the neutrino yield from W^W" pairs is similar to that for tt 
pairs. Therefore, neutrino rates are affected only weakly if the process xX ~^ WW* 
tWb is neglected. 

Annihilation in the Sun and Earth takes place only when the annihilating neutrali- 
nos have relative velocities v ~ 10~^ <^ 1. Therefore, if we expand the annihilation 



cross section as 



av = a + bv'^, (1.1) 

then we consider only the a (the s-wave) term. The rehc abundance also depends on 
the annihilation cross section. However, annihilation in the early Universe takes place 
when f ~ 0.5. In the early Universe, thermal averaging of the cross section smoothes 
out the jump in the annihilation rate near particle thresholds [^. Therefore, three-body 
final states should have only a negligible effect on the relic abundance, so we do not 
consider it further. 

We calculate the cross section and neutrino signal for the tt* final state in Section II 
and those for the WW* final state in Section III, and discuss these results in Section IV. 
Some lengthy matrix elements needed for the calculation are given in the Appendices. 



2. Calculation of the tWb final state 

We begin by calculating the cross section to the tWb final state. In some respects, our 
calculation parallels a recent calculation of the three-body decay of Higgs bosons into 
this final state The Feynman diagrams for this process are shown in Fig. |l]. Like 
annihilation to ti pairs, annihilation to this three-body final state takes place via s- 
channel exchange of and A° (pseudo-scalar Higgs) and t- and w-channel exchange of 
squarks. Although there are additional diagrams for this process, such as those shown 
in Fig. 1^, these are negligible for the regions of parameter space that we investigate for 
the reasons discussed in the Introduction. 




Figure 1: Feynman diagrams for annihilation of neutralinos to the tW b final state through 
an intermediate tt* state. 
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Figure 2: The contribution to the tWb final state from W^W pairs. These are negligible 
for the regions of parameter space we investigate. 



The total cross section for annihilation into three-body final states is given by 



2(27r)5Ai/2(s,rn2 m2)4 



llMl^d^PS; 



(2.1) 



where 



X{x,y,z)=x +y +z — 2xy — 2yz — 2zx ; 



(2.2) 



A4 is the matrix element for the process; Nc is the color factor; and d^PS is the 
three-body phase space which can be written as 
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P4, p5, and Pq are the four- momenta of the t quark, W boson and h quarks respectively, 
as shown in the Feynman diagrams; and a/s is the total center-of-mass energy. We 
approximate mf, = 0, and take the fj-ei limit where Wrei is the relative velocity of 



the annihilating neutrahnos. The boundaries of the phase space are given by 
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The total cross section is then given by 



1287r3 



/ dxe / 



^4max 



dx4-\M\ , 



where the spin-summed amplitude squared is given by: 

\M\' = J2\M\\ 



spin 



with 



(2.5) 
(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



M = Mit + M2t + M^t + Mu - Mi„ - M^u - M^^ - M^^ + M5 + M7. (2.11) 

If we make the approximation t,u ^ rrij (we cannot make this approximation for s in 
this case), then the MjS are given by 
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Here, the subscripts 4,5,6 denotes the t{i), W~(W~^) and b{b) respectively, and ps = 
Pb + Pe is the 4-momentum of the virtual top quark. The subscript t and u denotes t- 
and M-channel exchange of a sfermion. The left- and right-hand projectors L and R are 

iJ = l±25. L^l^; (2.12) 

and the coefficients c, are given by 

^ 2g^CriNi,-Ni,) ^ _ 2g\i{Nl^-Nl) _ ^ _ 2gTAoohAff . 
s cos^ 9w — 4:m'^ ' s cos^ 6'i„ — 4m^ ' s — 2m\ ' 

where the sum on j is over the squarks; j = 1, 2, 3 for left-handed squarks and j = 4, 5, 6 
for right-handed squarks. The quantities rrij, are the masses of sfermions j. The 
couplings Wj and are xf f right-hand and left-hand couplings respectively; they are 
defined as W^ ^ j q and X'u,t,j,o notation of fl]], and Taoo and /ia// are couplings 

of the Higgs boson A. N^q and N^q are the composition coefficients of the neutralino 
given in [|l|. Also, 

12 2 
ci = - sin^Ow; Cr = --sm'^9w (2.13) 



For the virtual-top propagator we use the relativistic Wigner-Weisskopf form 

^ + m 
p1 — rnl + irritTt 



(2.14) 



The decay width of the top quark is not yet measured, so we use the prediction of the 
standard model. 

We calculate using standard trace techniques. Notice also that 
Mu^jt = Miu,ju = Mij , M5,it = -M5,^u = M5,i , Mj^u = -My.i^ = Mr^i ; (2.16) 
where i,j = 1,2,3,4. Thus, we have 

\M\'^ = 2 [Mi,i + M2,2 + M3,3 + M4,4 + 2(Mi,2 + Mi,3 + Mi,4 + M2,3 + M2,4 + M3,4)] 
-2 [Mit,i, + M2t,2u + Mst,3u + M^^^u] + M5,5 + M7,7 
-4 [Mi<,2„ + Mit,3„ + Mit,4. + M2t,3« + M2t,4« + M3t,4„] 
+4 [M5,i + M5,2 + M5,3 + M5,4 + Mr,! + M7,2 + M7,3 + M7,4] . (2.17) 



These Mj ^ are given in Appendix |^. 



Ms = 260 - 400 GeV 


/i = 10^ GeV 


ruA = 10^ GeV 


tan/? = 2 


m? = lO^GeV^ 


m? = lO^GeV^ 



Table 1: Model parameters for Fig. y 
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Figure 3: The dashed curve is the total annihilation cross section times relative velocity to 
two-body final states. The solid curve includes the tWb final state as well. 



The three-body cross section for a series of typical models are shown in Fig. ^ 
The parameters of the models are given in Table |I|. Here we assume the GUT mass 
relation Mi = |M2 tan^ 9w, and the mass of the lightest neutralino varies as M2 varies. 
Note that our choice of the parameter 3> M2 ensures that the neutralino is primarily 
gaugino. We have also chosen diagonal and degenerate mass matrices for the squarks 
and sleptons. 

As expected, the three-body cross section approaches to the two-body value above 
the top threshold (we take rrit = 180 GeV). Below the top mass it is non-zero but drops 
quickly. 

We now consider the neutrino spectrum. If the rate for neutralino capture in the 
Sun is fixed, the flux of upward muon induced by energetic neutrinos from neutralino 
annihilation in the Sun and/or Earth depends on the various neutralino annihilation 
branching ratios, but not on the total annihilation cross section. The flux of upward 
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Figure 4: Detection rate of neutrinos from neutralino annihilation in the Earth versus 
neutralino mass. The ordinate is ^ defined in Eqs. ( |2.20| ) to which the detection rate is 
proportional. The dashed curve includes only two-body final states, the solid curve includes 
the tWh final state as well. 



muons (the quantity to be compared with experiment) is 

Tdetect OC ^total = ^ ^ Bf{Nz^) 



(2.18) 



where the sum on i is over muon neutrinos and anti-neutrinos; a, are neutrino scattering 
coefficients, = 6.8 and ap = 3.1; and 6j are muon-range coefficients, b,^ = 0.51 and 
bp = 0.67 [0. The sum on F is over all annihilation channels. The quantity z = Ey/rriy. 
is the neutrino energy scaled by the neutralino mass. In the two-body case, this is the 
same as the injection energy of the particle produced in neutralino annihilation, and 
the second moment is given by 



(2.19) 



The functions {N z'^) p ^{Ei^) have been calculated for injection of all particles that 
the neutralino may annihilate to in both the Sun and Earth for both muon neutrinos 
and anti-neutrinos |]^, ||. Note that the {Nz"^) are different for particles injected in the 



Sun and Earth since slowing of quarks and slowing and absorption of neutrinos may be 
significant in the Sun but not in the Earth. For the same reason, the {Nz'^) for muon 
neutrinos and anti-neutrinos will differ for particles injected into the Sun, although 
they are the same for particles injected into the Earth. 

In the three-body case, E^^ ^ m^. However, Eqs. ( p.l9| ) can be generalized. For 
the twh final state, the second moment is given by 

Bp{Nz') = ^ I dx,dx,\M\' {{Nz')^xl + {Nz')^xl + {Nz')^xl) . (2.20) 

We plot the result for annihilation in the Earth in Fig. ^. (The results for annihila- 
tion in the Sun are similar.) We plot the relevant combination of the second moments 
{Nz"^) given in Eq. ( ^.181 ). The predicted muon flux is proportional to this quantity. 
The dashed curves in the Figure show the results obtained by ignoring the tWb final 
state, and the solid curves show the results with the tWb final state. 

Although the three-body cross section is small except just below the ti threshold, 
its contribution to the the second moment Bp{Nz'^) may be important, as illustrated 
in Fig. ^. This is because the {Nz"^) for top quarks and W bosons is significantly 
larger than that for the light fermions P]. 



3. Neutralino annihilation to WW* 

Below the W^W~ threshold, the neutralino can annihilate to a real W and a virtual 
W which we denote by W*. The real and virtual W bosons then decay independently. 
About 10% of these decay into a muon (or anti-muon) and a muon anti-neutrino (or 
neutrino). These processes are illustrated in Fig. ^, where the W* decays into a fermion 
pair //', and the fermion pair can be rz/, yuz/, ev, cs, or ud. 

The WW* calculation is similar to the tWh calculation. We are only interested in 
the f — > limit. In this limit, only chargino exchange in the t and u channels shown in 
Fig. ^ are important. Neutrinos can be produced either by the virtual W* or by decay 
of the real W . 




Figure 5: Neutralino annihilation to Wff through a virtual W* by chargino exchange in 
the t oi u channel. 



The cross section is given again by Eq. ( p.9|) , and x^, x^, and Xq are defined similarly. 
Neglecting the fermion mass, the boundaries of the phase space can be written as: 

7714 



Xin 



X4r\ 



XQn 




where = mw- The amplitude is 

M = Mot + Mu - Mou - Ml, 



where 



i,t 



x7 



3fM 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 

(3.6) 
(3.7) 



Here i = 0,1 denotes the two chargino states, their masses, and Oi^i and Oi ji are 



the couplings. The quantity k = — p2 and k' = p^—pi are the momentum transfered 
in the t and u channel respectively, and in the f — limit, pi = p2 and k = k' . Also, 



9 



3/x 



y/2pl- + imwTw 



1 ) 



(3.8) 



and W^u is the wavefunction for the outgoing W boson. 
Straightforward calculation yields 



j=0 



{k^ - M^){k^ - Ml, 
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-E 

i=0 



(fc2 - M2)2 



where 



[{0t, + 0t^)M,^^ + 20^,0^^M,^,] 



9^{P5P6) 



{pI - rri^ f + m^r 



2 r2 ' 
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(3.9) 



-1 r\ 



Ma = 1000 GeV 


/i = 55 - 110 GeV 


TTiA = 1000 GeV 


tan /3 = 2 


m? = lO^OGeV^ 


mf = lO^OGeV^ 



Table 2: Model parameters for Fig. ^ 
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Figure 6: Cross section times relative velocity for neutralino annihilation to W jjLU. The 
dashed curve includes only two-body final states and the solid curve includes the three-body 
final states. 

and the matrix elements Mi^j are listed in Appendix 0. 

This cross section can be evaluated numerically. The result for models with param- 
eters given in Table |^ is shown in Fig. ^. Like the twh case, the three-body annihilation 
cross section drops quickly below niw and approaches zero as — > mw l'^- 

The calculation of neutrino flux is again analogous to the twh case. However, we 
will not consider the neutrinos produced by decay of the muon or other fermions as 
they are of negligibly low energy. We have 

Bp{Nz^) = J dx,dxe\M\' Q(l + ^/3^)x^ + x^^ /2, (3.10) 

where (3w is the velocity of the real W boson. The second term on the right-hand side 
accounts for the neutrino produced via the W*. The first term is from the decay of the 
real W boson. This W boson has a probability of Tw^f^u to decay to a muon neutrino. 
On the other hand this W boson is produced in all XX ~^ ^ff' channels, and each 
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Figure 7: Second moment of neutrino energy spectrum for neutrinos from neutralino 
anniliilation in the Earth. The dashed curve is the two-body calculation, and the solid curve 
includes the three-body final states from an intermediate WW* state as well. 



of these channels has approximately the same cross section, so the contribution has a 
weight of n^han^w^f^u ~ 1- 

The result for annihilation in the Earth is shown in Fig. 0. The results for annihi- 
lation in the Sun look similar. We see again the second moment agrees with the WW 
result above the threshold, but does not decrease immediately below mw, despite the 
fast drop of the cross section. In fact, the value of the second moment is slightly higher 
below the WW threshold, and does not drop until the cross section vanishes as 
approaches mw/^. Below mw/'^, the four-body channel might become significant, and 
may smooth the jump in much the same fashion as the three-body channel does near 
the two-body channel threshold, but we will not consider it here. 

To understand why there is a dip near the WW threshold in Fig. ^ consider Eq. 
( ^.10| ). There are two terms in the integrand. The term is the contribution from 
the decay of the "on-shell" W, and the Xq term is from the virtual W*. We plot these 
contributions separately in Fig. ^. 

To make comparison easier, what are actually plotted are twice the contributions 
from these terms plus contribution from other channels. As one can see from the plot, 
both the real and virtual W contributions agree with the two-body result above the 
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Figure 8: Contributions to the second moment ^ for annihilations in the Earth. Short 
dashed curve: two-body result; solid curve: three-body result; dotted curve: (twice of) real 
W contribution; long dashed curve: (twice of) virtual W* contribution. 

threshold as expected; below the threshold, however, their behaviors are quite different. 
While the virtual W* contribution smoothly drops, the real W contribution actually 
peaks well below the threshold; then it decreases as increases, and has a sharp turn 
at the WW threshold. This is because the the second moment is proportional to x1, 
and below the threshold steadily decreases as increases, as can be seen from 
Eqs. (p.l|)-(pT2D. However, above the threshold, is stabilized at the x^ = 1 pole, so 
it stops decreasing, and gradually increases as Pw increases. When averaged with the 
smooth curve of the W* contribution, this produces the strange dip at the threshold. 

In this calculation we have assumed that one W is on-shell while the other is off-shell. 
It could happen that both W^s are off-shell, and if such contributions are important 
this may produce a smoother curve and the dips in these plots may disappear. We 
examined this possibility by smearing the on-shell mass in the range mw — SF^y < 
1714 < mw + and calculating an average with the assumption that the cross section 
is proportional to 

1 

{ml - m^f - m\V%^' 

This produces a result that differs by less than 10~^ from the simple on-shell result. 
Hence it suggests that the on-shell approximation works well. 




4. Conclusions 



We made several approximations in our calculation. For example, in the tWb cal- 
culation, we neglected the contributions from diagrams involving production of the 
tWb final state through virtual W bosons. However, in those regions of parameter 
space where such omissions become unacceptable, the process we are considering is not 
very important. In the Wff' case, our calculation is generally applicable, but it is 
most useful for a higgsino, since in that case the W^W~ channel becomes large above 
threshold. 

Although the contributions of these three-body channels are important only in a 
limited region of parameter space, it may be a large effect. In fact, our calculation 
shows that although the cross sections of these annihilations are significant only just 
below the two-body channel threshold, due to the high energy of the neutrinos they 
produce, they can enhance the neutrino signal by many times and actually dominate 
the neutrino signal far below the two-body threshold. Furthermore, the regions in 
question may be of particular interest. For example, motivated by collider data, Kane 
and Wells proposed a light higgsino (but see []10|)- There are also arguments that 
the neutralino should be primarily gaugino with a mass somewhere below but near the 



top-quark mass |Tl 



There are many parameters in the minimal supersymmetric model. The results 
shown in Figs. § and ^ are of course model dependent, and these effects might be more 
or less important in models with different parameters. For a generic neutralino, these 
effects would be small. However, for the regions of parameter spaces we have focussed 
on it seems that the increase of the indirect-detection rate is fairly robust. For these 
regions of parameter space these effects should be considered for calculations of rates 
for energetic-neutrino detection. 
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A. The W cross section 

We give the Mj j for the tt* cross section calculation in this Appendix. 
For convenience, we introduce four-momenta 
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[these should not be confused with the cross-section factors introduced in Eq. ( |1.1|) ]. 
We further introduce 

1 (,2 

We can then write the Mj^, as 

Mi,i = mclm^t{pib){p2P4) , M2,2 = ^Nclpl{pia){p2P4) , 
M3,3 = ANclml{pih){p2PA) , M4,4 = ^Nclpl{pia){p2Pi) , 



j I I 

M5,5 = 16A^ {C5P3 [{piPi){p2a) + (pia)(P2P4) - "^^(P^a)] H 
+C6m3 [(piP4)(p2&) + {pi^){p2'Pi) - rn^ipih)] 

+2c5C6m2(p36) - 2mJ)) } , O 

M7,7 = 16iV c?m2 {2 [(P3P4) + (P36) + [ml - p^] (p^fo) } . ^ 

The Mi/s are given by I ^ 

Mi,2 = ANciC2mlmf{psb) , Mi,3 = 4A^ciC3m^m^(pi6) , <X> 

Mi,4 = ANciCAm^mt{p2h){p2Pi) , M2,3 = '^Nc2C^m^mt{pzh){p2Pi) , ^ 

M2,4 = ANc2Cim^mt{p^p:i){pia) , M3,4 = 4:Nc3C4mlmf{p3b) , 00 

and ^ 

Ms^i = -ANciml {c5 [(piP2)(P3&) + {.P2Pz){.Pib) - {,Pm)ip2h) - 2mJ(p36)] Q 

+C6 [2(^2^4) - } , I — ^ 



M5,2 =4A^C2 {C5 [4(^2^4) (PlP3)(P3&) " 2(^2^4) (PsPs) 

-2mJ(p3P4)(p36) + mJ(p3P3)(p46)] 

+C6m^ [(piP2)(P3&) + iPiP3){P2b) - {p2Pz){pih) - 2m\{p^h)\ ] , 

M5,3 = iNcsmtm^ {2C5 [(P2P4)(P3&) + {PlPi){P2.b) + {PlP3){Pib) - {p3P4){pib)] 

\c^m\ [(p3&) -2(pi6)]}, 

^5,4 = - 4A^C4m;^ {C5(P3P3) [2(^2^4) - (P1P4)] + Cf,m\\{p2p-i) - 2(piP3)]} , 

My^i = 2NciC7m^mt (2(^36) [(^3^4) +772^2] + {p4b)[ml - pj]) , 

M7,2 = 8Nc2C7mlmt{p3b), Mj^s = AN c^C7mlm\{pih) , 

M7,4 = 2A^C4C7 {pl{piP2){PAa) + ml{pip2){p?.h) + mlm^ipsb)) , 

M7,5 = 8A^C7 (4c5mimJ(p36) - Cgmtm^ [2 ((^3^4) + mDip^b) + {ml - pl){pib)]) . 



Mit,2u = 2NciC2m^ml{psb) , 



M2t,2u = 2Nclmlpl{p4a) , 



Mu,iu = 2Nclmlml{pih) , 
Mit,3u = 2NciCzmlm^{p2h) , 
M2t,Au = 2Nc2Cimtm-xpl{p2a) , M3t,4.u = 2A/"c3C4mt (pip2)(p3&) , 

Mit,4u = '^NCiC4mtm-^[{p2P4){P3b) + (P2P3)(P4&) - (P3P4)(P2&)] , 
M2t,3„ = 2Nc2C3mtm^[{p2P4){p3b) + (P3P4)(P2&) - (P2P3)(P4&)] , 
M3t,3u = 2A^C3m?[(p2P4)(pi&) + (P1P4)(P2&) - (piP2)(p4&)] , 

-/V^4t,4« = '^Nclpl[{p2P4){pia) + {piP4){P2a) - {piP2){P4a)] . 



In the f — > hmit, we can approximate pi = P2 = 0. In terms of X4 and Xe, we 
have 
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X4 + 


1 . , 

^^645) , 
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1 . , 
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— X4 — -A564j , 



where 



^i,j,k — 



mf + mi 
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mt 
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B. The WW* cross section 

We give the M^j necessary for the WW* calculation in this Appendix. 
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As in the tWb case, these four-vector products can be expressed by the variables 
X4, X5 and Xq: 



{P1P2) ^ m\ , 

(piPe) ~ {P2P&) ~ m^xe, 
(P4P5) ~ 2m^(l - xq 



(piPs) ~ (P2P3) ~ m^(2 - X4) 
(piPs) ~ (P2P5) ~ m^xs , 
P3 ~ 2(p5P6), 



m. 



m^ 



4m^' 



(P4P6) ~ 2m (1 - X5 + — ^) 



4m^ 



Ami' 



and note 



k = Pi- p2, 

SP5Pl)Pe + iP6Pl)P5 



(PbPe) 



(B.l) 
(B.2) 
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